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On  Nonlocal  Continuum  Mechanics 


by 

A.  E.  Green^  and  P.  M.  Naghdi* 


Abstract.  It  is  shown  in  this  paper  that  the  conservation  laws  of  nonlocal 
continuum  mechanics  are  equivalent  to,  and  can  be  derived  from,  the  conserva- 
tion of  energy  and  the  invariance  conditions  under  superposed  rigid  body 
motions.  Also,  the  theory  of  nonlocal  thermoelasticity  is  reconsidered  in 
the  light  of  recent  developments  in  thermodynamics,  taking  invariance 
conditions  fully  into  account. 


1. 


Introduction 


In  the  mechanics  of  single  phase  continue  in  which  nonlocal  effects  are 
excluded,  it  is  known  that  the  equations  of  conservation  of  mass,  momentum  and 
moment  of  momentum  may  be  derived  from  the  energy  equation  with  the  help  of 
invariance  conditions  under  superposed  rigid  body  motions.  In  the  first  part 
of  the  paper,  we  prove  that  the  same  is  true  for  nonlocal  continuum  mechanics. 

From  another  point  of  view,  this  shows  that  there  is  consistency  between  the 
equations  of  conservation  of  mass,  momentum  and  moment  of  momentum,  and  the 
equation  of  conservation  of  energy  if  one  accepts  the  invariance  criteria. 

Constitutive  equations  for  a nonlocal  elastic  and  thermoelastic  solid 
have  been  studied  by  a number  of  writers  (see,  e.g.  [4,6,7]  and  references 
contained  therein).  There  appears  to  be  a number  of  differences  between  these 
authors  over  invariance  conditions  under  superposed  rigid  body  motions  and 
subsidiary  conditions  to  be  satisfied  by  force,  couple  and  energy  residuals. 

Recently  an  entirely  new  approach  to  continuum  thermodynamics 
has  been  proposed  by  Green  and  Naghdi  [8].  In  view  of  the  differences  between 
previous  authors  and  in  the  light  of  the  developments  in  [8],  constitutive 
equations  for  a nonlocal  thermoelastic  solid  are  reconsidered  here  taking  the 
invariance  conditions  fully  into  account. 

Consider  a finite  body  B with  material  points  X and  identify  the  material 
point  X with  its  position  X in  a fixed  reference  configuration.  A motion  of 
the  body  is  defined  by  a sufficiently  smooth  vector  function  x which  assigns 
position  x = x(X,t)  to  each  material  point  X at  each  instant  of  time  t.  In  the 

mm 

present  configuration  at  time  t,  the  body  (B  occupies  a region  of  space  R bounded 
by  a closed  surface  3R.  Similarly,  in  the  present  configuration,  an  arbitrary 
material  volume  of  <B  occupies  a portion  of  the  region  of  space  R,  which  we 
denote  by  a part  p (eft)  bounded  by  a closed  surface  dp.  Let  p=  p(X,t)  be  the 
mass  density  in  the  present  configuration  and  designate  the  velocity  vector 

-i 


2. 


at  tine  t by  v-x,  where  a superposed  dot  stands  for  the  material  time  derivative. 

/v 

Throughout  the  paper,  we  assume  that  mass  is  conserved  locally. 

We  assume  that  in  its  present  configuration,  the  body  ® is  acted  upon  by 
an  external  body  force  b = b(X,t)  per  unit  mass  and  an  external  surface  force  t 
per  unit  area  acting  on  the  boundary  dR.  The  rate  of  work  of  these  external 
forces  are,  respectively,  b • v per  unit  mass  and  t . v per  unit  area.  We  also 

^ CW  A/ 

admit  the  existence  of  an  internal  body  force  f=f(X,t)  and  an  internal  body 
couple  represented  by  an  axial  vector  \=\(X,t),  each  per  unit  mass,  with 
corresponding  rates  of  work  per  unit  mass  equal  to  f • v and  §-  X • w,  respectively, 
where  w=curl  v is  the  (axial)  vorticity  vector.  In  addition,  there  is  an 

^ /v 

internal  surface  force  t = t(X,t;n),  acting  over  each  surface  dp  with  outward 
unit  normal  n,  such  that  t • v represents  rate  of  work  per  unit  area;  the  field 
t,  called  the  stress  vector,  is  measured  per  unit  area  of  dP  and  assumes  the 
value  t on  dR. 

t 

WiVth  reference  to  the  thermal  properties  of  the  body,  we  first  introduce 
the  absolute  temperature  at  each  material  point  by  a scalar  field  0=  9(X,t)>0, 


as  well 


s an  external  rate  of  supply  of  heat  r = r(X,t)  per  unit  mass  and  an 


external  rate  of  surface  supply  of  heat  -h  per  unit  area  acting  across  dR. 

Also,  we  admit  the  existence  of  an  internal  surface  flux  of  heat  -h  =-h(X,t  ; n) 

m /v 

per  unit  area  across  the  surface  dP;  the  field  h,  called  the  heat  flux  and 
measured  per  unit  area  per  unit  time,  assumes  the  value  h on  dR.  We  define 
the  ratios  of  r and  h to  temperature  by 

s = s (X  ,t ) = “ , k = k(X,t;n)  = £ (l.l) 

and  refer  to  (l.l).  0 as  the  external  rate  of  supply  of  entropy  per  unit  mass 


and  the  internal  surface  flux  of  entropy  per  unit  area  of  dP,  respectively; 
the  latter  assumes  the  value  k = h/0  on  the  external  bouruary  of  ffl,  where  k is 


3 


the  external  rate  of  surface  supply  of  entropy  per  unit  area  of  dR.  Further, 
we  assume  the  existence  of  the  specific  entropy  specified  by  a scalar  field 
T)=  *n(X,t)  and  an  internal  rate  of  production  of  entropy  £=§(X,t)  per  unit  mass 

A# 

with  a contribution  to  the  internal  rate  of  production  of  heat  equal  to  0§  per 
unit  mass.  In  addition,  we  admit  the  existence  of  a potential  e=e(X,t)  called 
the  specific  internal  energy  and  a residual  rate  of  production  of  energy 
0=#(X,t)  per  unit  mass. 


2 . The  energy  equation.  Invariance  conditions. 

The  balance  of  energy  for  an  arbitrary  material  volume  occupying  a part 
in  the  present  configuration  can  be  stated  as 


— - l p(e  + | v • v )dv 

& L Wp  r+J 


= f p(r  + b-  v + f*v  + i\  *w  + 0)dv+  I (t.v-h)da 

JP  ~ Jc>P 


Under  superposed  rigid  body  motions  represented  by 


= c +Qx 


where  c(t)  is  a vector  function  and  Q(t)  a proper  orthogonal  tensor  function 
of  time,  we  assume  that  the  mass  density  p,  the  heat  supply  r,  and  the  fields 


T) , 9,€  ,0  ,h ,f , X , t , 

rv 

nl 1 are  unaltered  apart  from  orientation  in  the  case  of  vectors.  Thus,  if 
T\,9,P,r,e  ,. . . ,t  become  T)4 ,9+, . . . ,t+  under  superposed  rigid  body  motions,  we 
have 

T)+  = T|  , 0+  = e , 

+ + + + , + , 
p=p  , r = r , e =*  c , 0=0  > h=h  , 

f+  = Qf  , X+  = Q X , t+  = Qt  . 

In  addition  to  (2.1),  which  holds  for  every  part  P (eft),  we  assume  that  the 
energy  equation  for  the  whole  body  is  the  same  as  that  for  local  continuum 


P 


(2.1) 


(2.2) 


(2.3) 


(2.4) 


mechanics  so  that 


I (e  + 5 v • v)dv  = f p(r  + b • v)dv  + (*  (t  • v - h)da  . 
Jft  ~ ~ ~ ~ JdR  ~ ~ 

It  follows  from  (2.1)  and  (2.5)  that 


(2.5) 


I p(f  • v + \ \ • w + 0)dv  = 0 . 

JR 

For  local  continuum  mechanics,  a derivation  of  field  equations  and  jump 


(2.6) 


conditions  from  the  energy  equation  has  been  given  by  Green  and  Rivlin  [l]  and 
by  Green  and  Naghdi  [2].  Here  we  employ  a method  of  derivation  wnich  is  modelled 
on  that  used  by  Naghdi  [3,  p.  487]  in  a different  context.  Let  Cl  be  a second 

r-Vy 

. rp 

order  tensor-valued  function  of  time  defined  by  Cl(t ) = Q(t)Q(t)  , where  the 
T 

notation  Q stands  for  the  transpose  of  Q.  Then,  since  Q is  a proper  orthogonal 
tensor,  we  have 


Q = o Q , ci  = - ci 


(2.7) 


Hence,  there  exists  an  axial  vector  uj  corresponding  to  the  skew  tensor  Cl  such 


Cl  a = a)  x a 


(2.8) 


for  all  vectors  a.  Naghdi  [3,  p.  453]  has  shown  that  Q can  be  chosen  so  that 


n = n , u>  = , 


(2.9a) 


s.<t) » i . « * a,  . 


(2.96) 


where  Cl  is  a constant  skew-symmetric  tensor,  ou  is  a constant  vector  and  I is 

r*JQ  *** 

the  identity  tensor.  With  this  choice,  in  all  subsequent  work  we  may  write 


0 u>  = 0 , Qu)  = 0 , Cl  a . a = 0 


(2.10) 


for  all  vectors  a . 


6. 


Consider  first  a motion  of  the  continuum  which  differs  from  the  given 
motion  only  by  a uniform  rigid  body  translational  velocity  specified  by 


x = x + c , c = 0 , 


+ • + 
v = v + c , w = w 


From  (2.4)  and  (2.11)  it  follows  that  p,r ,e ,0,h,f and  t are  all  unchanged 

under  this  superposed  rigid  body  motion.  We  assume  that  the  new  motion  can  be 

maintained  by  the  same  body  force  b.  The  energy  equation  associated  with  this 

• • 

motion  is  obtained  from  (2.1)  by  replacing  v by  v + c,  where  c is  a constant 

esj  /w  /■w 

vector.  Then,  remembering  that  the  volume  element  dv  and  surface  element  da 
are  also  unchanged,  by  subtraction  we  obtain 

f -jx  I p dv  + c • [jj-  f pv  dv  - [ p(b  + f )dv  - [ t da]  = 0 

~dtJp  ~ dtjp~  ~ ~ ~ 

for  all  arbitrary  c,  the  quantities  in  the  brackets  being  independent  of  c. 
Hence , 

AJ^pdv.o  , (a 

f pv  dv  = f p(b  + f)dv  + f t da  (2 

dt  Jp  ~ ~ ~ Jap  ~ 

These  are  the  equations  of  conservation  of  mass  and  linear  momentum.  From 
(2.12)  and  (2.13),  under  suitable  smoothness  assumptions , follow  the  field 
equations 


p + p div  v = 0 


t = T n , div  T + p(b  + f)  = pv  , 


where  T is  the  stress  tensor. 


7. 


By  applying  the  same  invariance  conditions  to  the  energy  equation  (2.5) 
for  the  whole  body,  we  obtain 


d (*  . r r 

— pv  dv  = | pb  dv  + t da 

JR  ~ JR  ~ JdR  ~ 


JR  ~ JR  °dR 

From  comparison  of  (2.15)  with  (2.13),  we  conclude  that 


(2.15) 


i pf  dv  = 0 
JR  ~ 


(2.16) 


The  condition  (2.16)  also  follows  from  (2.6)  by  a similar  consideration  of 
invariance.  It  states  that  the  total  internal  force  acting  on  the  whole  body 
is  zero. 

Next,  we  consider  a superposed  rigid  body  motion  specified  by 


x = Q x 


(2.17) 


and  restrict  our  attention  to  the  case  in  which  the  conditions  (2.9)  hold. 


8 

Then, 


v = Q v + Q Q x , 


•+  • 2 + 
v = v + 2Uv  + Qx  , w =w+2o)  , 

*>•*  rs^  rZ /w  ~ Q 


(2.18) 


and 


+ + 


+ + 


f - v = f.v+u)_-xxf  , X.  * w = X • w + 2ou  . X.  , 

^ rv#  /v  rvQ  /v  (V  sns  rw  rv  isO  /vi 


t+  • v+  = t • v + uu  • x x t , 

/v  /v  /Vi  /vO  IV  /V/ 


(2.19) 


+ + T 2 

v • v = v • v + 2Q  id  ■ xxv-fl  • Q x Q x , 

i-v  v iv  /v  /v  rvQ  iv  /v  <vO  /v  xv  rw  iv 


where  in  obtaining  (2.19)  use  has  been  made  of  the  requirements  (2.4).  We  now 
need  to  consider  the  invariance  condition  to  be  satisfied  by  the  body  force  b 


^The  result  for  v+  is  recorded  without  the  use  of  (2.9b)  since  v+  • v+  occurs 

^ iv  iv 

later  in  an  integral  which  is  differentiated  with  respect  to  time. 


1 


under  superposed  rigid  body  motion  (2.17),  i.e.,  we  must  find  the  new  body 
force  b which  will  sustain  the  motion  (2.17).  From  (2.14)  and  (2.4),  we 
obtain  the  relation+ 


8. 


b+  - v+  = Q(b  - v) 


(2.20) 


and  the  rate  of  work  of  body  force  b+  per  unit  mass  due  to  v+  is 

b . v = b • v + u)  •xxb-in*(vx  + xv)  , (2.21 ) 

where  (2.9) > (2. id)  and  (2.2 0)  have  been  used. 

Consider  now  an  energy  equation  of  the  form  (2.1)  associated  with  the 
motion  (2.17)  and  subtract  from  it  the  energy  equation  (2.1).  Then,  with 
the  help  of  (2.10),  (2.19)  and  (2.21),  it  follows  that 

ou  • [~-  I px  x v dv  - I p(x  x (b+f ) + \)dv  - i X X t da] 

~0  dtJp~~  Jp  ~ J^p  ~ ~ 

- 2 Cp  • [-^r  [*  p(Q  x)(Q  x)dv  - f p(v  x+x  v)dv]  = 0 

r>UJ  CIO  J ^ /v  *■*-»  J ^ A#  a/  rsJ  r+s 

for  all  arbitrary  ou^.  From  the  last  result  follow  the  two  equations 

-7T  f pxxvdv  = f p[x  x (b+f ) + A.]dv  + [ xxtda  , (2.22) 

dt  ~ ~ Jp  ~ — ~ JaP  ~ ~ 

TT  f P(Q  x)(Q  x)dv  = f p(v  x + x v )dv  , (2.23) 

the  first  of  which  represent  the  conservation  of  moment  of  momentum.  With  the 
use  of  (2.14),  the  field  equation  resulting  from  (2.22)  can  be  reduced 
to 


pr  + T - T = 0 , 

r r+s  A/  A^  /SJ 


(2.24) 


Alternatively,  the  result  (2.20)  can  be  obtained  directly  from  (2.13)  and 
(2.14);  in  this  connection,  see  [ 3 > P*  486]. 





9. 


where  T is  the  skew  tensor  corresponding  to  the  axial  vector  \ so  that 

£ a = Xxa  (2.25) 

for  all  arbitrary  vectors  a. 

Suppose  now  that  there  are  no  surfaces  of  discontinuity  in  p.  Then,  in 
view  of  (2.9)  and  (2.14)^,  the  local  equation  resulting  from  (2.23)  is 
satisfied  identically  and  yields  no  new  information.  The  same  conclusion 
holds  even  if  the  mass  density  p has  a jump  discontinuity  across  a surface  P, 
assuming  that  x is  continuous. 

If  we  consider  th  * superposed  rigid  body  motion  (2.1?)  and  use  the  energy 
equation  (2.5) > by  a procedure  similar  to  that  which  led  to  (2.22)  we  find 

dr  r p 

— I px  x v dv  = ] px  x b dv  + I x x t da  . (2.26) 

dt  JR  ~ ~ Jft  ~ ~ Jdft  ~ ~ 

Comparison  of  (2.26)  and  (2.22)  implies  that 


(• 


p(x  x f + X)dv 

r+j 


0 . 


(2.27) 


The  condition  (2.27)  may  also  be  deduced  from  a similar  consideration  of 
invariance  of  (2.6)  under  the  superposed  rigid  body  motion  (2.17).  It  states 
that  the  total  internal  couple  acting  on  the  whole  body  is  zero.  We  observe 
that  (2.16)  can  also  be  derived  from  (2.27)  with  the  help  of  invariance 
conditions.  If  we  make  use  of  (2.24)  and  (2.25)  the  condition  (2.27)  can  be 
put  in  the  alternative  form 

[ [ p{x®f  - f 8>x)  + T- TT]dv  = 0 . 

JR 


(2.28) 


For  later  convenience  we  set 


|A  = 0 + f.V  + |-X-W 


L 1 


(2.29) 
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so  that,  from  (2.6),  we  have 

J pp.  dv  = 0 . (2.30) 


In  view  of  the  invariance  conditions  (2.4),  (2.1l)  and  (2.18),  it  follows  that 
H is  not  unaltered  by  superposed  rigid  body  motions.  The  energy  equation  (2.1) 
in  which  we  use  the  substitution  (2.29)  for  together  with  (2.5)  and  the 
derived  equations  (2.16),  (2.12),  (2.13),  (2.15),  (2.16),  (2.22),  (2.26), 
(2.27),  (2.28)  and  (2.30)  are  equivalent  to  those  postulated  by  Edelen,  Green 
and  Laws  [4]  and  Green  and  Naghdi  [5],  apart  from  notational  differences. 

They  are  also  similar  to  those  employed  by  Eringen  and  Edelen  [6]  and  Eringen 
[7j  except  that  these  authors  allowed  for  nonlocal  mass  changes.  However,  m 
[6],  a different  view  from  that  given  here  is  adopted  for  the  invariance 
condition  to  be  satisfied  by  the  function  which  corresponds  to  u,. 


11. 


3*  Balance  of  entropy 

Following  the  recent  approach  of  Green  ana  Naghdi  [8]  to  thermodynamics, 
we  assume  a balance  of  entropy  in  the  form 


d.  i*  P P 

— i pT|  dv  = | p(s+|)dv  - k da 

JP  JP  JdP 


(3-D 


for  every  material  volume  occupying  a part  p in  the  present  configuration, 
where  s and  k are  defined  by  (1.1 ).  The  field  equation  resulting  from  (3.1) 
is 


pT)  = p(s+5)  - div  p , (3.2) 

where  the  entropy  flux  and  the  heat  flux  vectors  p and  q are  related  by 

q=0p,k  = p.n,h  = q-  n.  (3*3) 

Invariance  conditions  additional  to  those  in  (2.4)  are 

k+  = k , p+  = Q p , q+  = Q q , 5+  = ? , s+  = s . (3.4) 

With  the  help  of  (2.14),  (2.29) > (3-2)  and  (3«3)>  the  field  equation 

resulting  from  the  energy  equation  (2.1)  can  be  reduced  to 

- p(lf  + Vfi)  + T • grad  v - pf  • v + pp 

- P?0  - P • grad  0 = 0 , (3.5) 

where  the  specific  Helmholtz  free  energy  is  defined  by  j/  = c - T]0. 

In  view  of  (2.4)  and  (3*4)  it  is  clear  that  p cannot,  in  general,  be  invariant 
under  superposed  rigid  body  motions  if  the  energy  equation  is  to  remain 
invariant.  The  behavior  of  p under  such  motions,  which  is  determined  by 
(2.4),  (2.24)  and  (2.29),  is  precisely  that  which  will  render  (3-5)  invariant. 


12. 


Once  constitutive  equations  have  been  specified  for 

ii  »£»j4,p»T  (3-6) 

the  energy  equation  (3-5)  is  regarded  as  an  identity  for  all  processes.  We 
also  regard  (2.24)  and  the  subsidiary  conditions  (2.16),  (2.28)  and  (2.30)  as 
identities  to  be  satisfied  for  all  thermo-mechanical  processes,  i.e.,  for  all 
motions  and  all  temperature  distributions.  As  far  as  the  last  four  of  these 
subsidiary  conditions  are  concerned,  the  view  taken  here  appears  to  differ 
from  that  adopted  in  [6,7]. 

When  the  body  is  in  equilibrium  with  v = 0 and  all  constitutive  variables 
independent  of  the  time  (but  dependent  on  X),  we  assume  that  0,  and  hence  n, 
vanishes.  Then,  it  follows  from  (3*5)  and  (3-2)  that 


p£9  + P * grad  0=0  , 

(3-7) 

p(s+§)  - div  p = 0 , 

pr  - div  q = 0 , (3-8) 


It  was  observed  in  [8]  that  the  quantities  £,Tl,i|f  may  be  arbitrary  to 
the  extent  of  additive  functions  f,f,-0f,  respectively,  where  f is  an  arbitrary 
function  of  x>9,  their  space  and  time  derivatives  and  a functional  of  their 
past  histories  throughout  R.  The  additive  functions  have  the  property  that 
they  make  no  contribution  to  the  differential  equations  for  Xj9  or  the  energy 

/v 

identity  (3*5 )>  boundary  and  initial  conditions.  The  arbitrariness  was  removed 
by  setting 

f = 0 . (3.9) 

The  same  situation  holds  for  the  present  nonlocal  theory  except  that  now  f may 
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also  be  a functional  of  X,0  over  the  whole  region  R.  The  arbitrariness  is  again 
removed  by  the  specification  (3. 9)* 

A mathematical  intepretation  of  the  second  law  of  thermodynamics  was  shown 

in  [8]  to  lead  to  an  inequality  which  reflects  the  fact  that,  for  every  process 

associated  with  a dissipative  material,  a part  of  the  external  mechanical  work 

supplied  to  a material  volume  P is  always  converted  into  heat  and  cannot  be 

withdrawn  from  p as  mechanical  work.  We  assume  that  the  constitutive  response 

functions  for  e,T]  include  dependence  on  the  set  of  variables  v',0',^'  and  their 

higher  space  and  time  derivatives,  where  a superposed  prime  denotes  the  values 

of  the  variables  at  any  point  in  the  present  configuration  R [see  Eq.  (4.21)]. 

* * 

We  refer  to  this  set  of  variables  collectively  as  Uj.  Further,  let  c ,T]  denote 
the  respective  values  of  c ,1]  when  the  variables  le  are  set  equal  to  zero  in  the 
response  functions.  Thus,  for  example, 

c = *(F/,x  ,0  ,Uj)  , « =«(F  ,x  ,0  ,0)  , U>=  (v  ,0  ,g  , . . . ) , (3*10 ) 

rsj  /w  M M 

where  the  dots  in  (3.10)^  refer  to  the  higher  space  and  time  derivatives  of 
v' , 0',g'.  Then,  the  inequality  for  nonlocal  continuum  mechanics  corresponding 
to  that  obtained  previously  for  local  continuum  mechanics  (see  (4.3)  of  [8])  is 


where 


[ pwdv  i 0 , 

R 


(3-11) 


pw  = - p(i|  + T\*Q)  + T • grad  v - pf  v = p[f  - j ♦ (‘H-T]  )o]  + P5®  + P • g-  PH  » 

* # * 

♦ = « - 07)  • 


(3.12) 


In  addition  to  a thermodynamic  restriction  corresponding  to  (3. 11),  two 
further  inequalities  which  place  restrictions  on  the  heat  conduction  vector 
and  internal  energy  are  proposed  in  [8,  Sec.  4].  But  the  corresponding 
additional  inequalities  will  not  be  needed  in  the  present  paper. 
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4.  Nonlocal  thermoelasticity 

As  noted  in  section  1,  a number  of  writers  have  studied  nonlocal  elasticity 
and  thermoelasticity  in  the  context  of  nonlinear  continuum  mechanics.  Without 
giving  an  extensive  bibliography,  we  refer  to  [4,6,7]  and  references  cited 
therein.  An  examination  of  these  papers  shows  that  there  is  no  clear  agree- 
ment on  the  constitutive  theory  of  thermoelasticity.  By  way  of  additional 
background,  we  recall  that  a fairly  general  class  of  thermoelastic  solids  was 
considered  in  [4];  and,  although  this  work  was  not  the  most  general  possible 
within  the  framework  of  the  nonlocal  theory,  the  results  given  satisfy  all  the 
invariance  conditions  and  subsidiary  equations  in  a manner  which  is  in  line 
with  those  of  the  present  paper.  More  general  constitutive  assumptions  were 
made  initially  in  [6,7],  but  the  results  in  [6,7]  in  the  case  of  zero  heat 
conduction  appear  to  be  different  from  each  other,  partly  because  of  the 
difference  of  interpretation  of  invariance  conditions. 

We  introduce  further  notations  by 

• 

g = grad  9 , F = Grad  X > F = L F , L = grad  v , 

zv  'W  rw  rw  rw  /v 

(4.1) 

I=(p/p0}££T  ’ T.  L = (p/po)P.F  , g = grad^-LTg  , 

where  the  operator  grad  is  with  respect  to  x,  the  operator  Grad  with  respect 
to  X and  p is  the  density  of  the  solid  in  a reference  configuration.  Also, 

p=  p(x,t)  , e=e(x,t)  , g = g(x,t)  , x = x(x,t)  , F = F(x,t)  , v = x(x,t)  , 

rw  «V  /w  M r+j  M M /w  o*  r>j  A#  ^ 

p'=  p(x',t) , e'  = e(x',t)  ,g/  = g(x',t)  ,x'  = x(x',t) , F'  = F(x',t) , v'  = i(x',t)  , 

(4.2) 

P0  = po(£}  ’ po  = Po(£'}  ’ dV  = dV(£}  » dV'  = dV<*'>  > 

dv  = dv(X,t)  , dv'  = dv(X',t) 


- ■ 
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In  conformity  with  (4.2),  henceforth  we  use  a superposed  prime  to  designate 
the  operation  of  replacing  X by  X'  and  X'  by  X.  For  example,  corresponding  to 

f = fCQjgjXjF.e'.g'.x'jF'.XjX')  , 

A#  M /V  M M M 


we  write 


f'  = f(0/,g',x#,F',0,g,x,F,X',X) 


For  a nonlocal  thermoelastic  solid,  we  adept  the  constitutive  assumptions 


are  functions  of 


9>£>*»£>£ 


and  functionals  of 


e'jg'jx'.F'jX' 


over  the  whole  region  R.  In  equilibrium,  the  deformation  function  x and  the 
temperature  0 on  the  right-hand  sides  of  ( 4.2 ^ and  (4.2)g  are  independent 
of  t but  can  be  chosen  arbitrarily.  The  function  p,  vanishes  in  equilibrium 
and  hence  vanishes  when  the  functions  in  (4.4)  and  (4.5)  depend  also  on  t. 
Subsequent  investigations  shows  that  p = 0 yields  an  inadequate  description 
for  nonlocal  thermoelasticity.  Inspection  of  the  energy  equation  (3-5)  as  an 
identity,  particularly  the  term  \|r,  suggests  that  we  allow  p to  be  a linear 
function  of  degree  one  in 


9»g,F,v 


(4.6) 


.... 
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and  a linear  functional  of  degree  one  in 


e\g',F>' 


over  ft,  as  well  as  depending  nonlinearly  on  (4.4)  and  (4.5). 

Before  proceeding  further  and  in  order  to  simplify  the  rest  of  our 
discussion,  it  is  convenient  to  dispose  of  some  mathematical  preliminaries. 
For  this  purpose,  we  first  introduce  the  vector  space  \j  of  all  quadruplets 
(a,a,b,A)  consisting  of  the  scalar  a,  the  two  vectors  a,b  and  the  tensor  A. 
The  space  V is  a 16-dimensional  vector  space  with  the  inner  product  of  the 
two  quadruplets  (a,a,b,A)  and  (c,c,d,B)  as 

r>s  ****  rss 

(a,a ,b,A)  < (c,c,d,B)  = ac  + a-c  + b-d  + A-B 

M M M M M M M M M M 

and  with  the  norm 

i 

| (a  ,a >b,A) I = [(a,ajb,A)  • (a5ajbjA)]^  • 

Let  U be  an  open  and  connected  subset  of  1/  consisting  of  those  quadruplets 
(a,a,b,A)  with  a>0  and  det  A>0.  Then,  ( 9,g,X >F)  certainly  lies  in  U.  We 

M M M M M IN/  M 

assume  all  variables  in  (4.4)  to  (4.7)  are  continuous  functions  of  X in  ft 

M 

for  all  times  and  introduce  the  abbreviation 

A(  • ) = (e(-),g(-),x(- ) >F( • ) ) • 

As  will  become  apparent  shortly,  the  constitutive  relations  to  be  introduced 
below  lead  us  to  the  study  of  the  functional  relations  of  the  form 


f(X)  = f ( A(X ' ) , A(X)  ) = f(A(*),A(X))  , 

M MM  MM  'V  M M 

X'  6ft. 


(4.7) 


(4.8) 


(4.9) 


(4.10) 


(4.11) 


where  dependence  on  t is  understood  and  where  ftQ  is  the  region  occupied  by  the 
finite  body  ffl  in  its  reference  configuration. 
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Let  v(*)  be  a given  influence  function  which  is  positive,  monotone 


decreasing  and  continuous  in  and  introduce  the  norm 

||r(*)||  = [f  y(|x-x'| )r(x') . r(x,)dv']^  (i 

Ro 

for  all  continuous  functions  T(*)  defined  for  every  X'£R  whose  values  lie  in 

~ ~ o 

V and  for  definiteness  also  specify  y(0) = 1.  The  collection  of  all  such  func- 
tions £(•)  forms  a pre-Hilbert  space  in  which  the  inner  product  for  the  two 
functions  r,  ( • ) and  ' • ) is 

<rL(-),re(-)>  = J v( |x-x ' | ' ) • r^(x')dv'  . (4 

R 

o 

We  suppose  that  the  functional  f(r(*)»A)  is  defined  for  every  A in  U and 
for  every  function  T( • ) in  H.  Then,  f is  said  to  be  continuous  if  for  each 
fixed  r(-)  and  for  every  y(«)  in  H, 

f(r(*)  + v(-),  A + n)  - f(r(-),A)  - o (4 

M /w  ev  /V 

as  |rr|  + ||v(‘)||  -•  0.  The  functional  f is  continuously  differentiable  if  the 
relation  (4.l4)  can  be  replaced  by 


MTTR01I  f <£<  • ) * v(  • ) , A . n)  - f (£(  - ) ,A) 


- Df(r(-),A)  • n-  6f(r(-),A|v(-))  - 0 

^ ^ /V  (V 1 /w 

as  I tt|  + ||v(*)||  -*  0.  Here  Df(r(*)»A)  is  a continuous  functional  with  values 
in  the  space  U and  is  simply  the  partial  derivative  of  f with  respect  to  A 
holding  r(-)  fixed,  i.e., 


Df(r(*),A)  = I7  (r(-),A)  . (4. 

~ ~ oA  ~ ~ 

/w 

Also,  for  each  fixed  r(*)  and  A,  6f (T( • ) >A|  v(  • ) ) is  a continuous  scalar-valued 
functional  defined  in  M which  is  linearly  dependent  on  v(')  an(i  is  continuously 


-t*: 
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dependent  on  the  pairs  (r(*)»A).  If  the  functional  f is  continuously  different!- 
able  in  the  sense  just  described,  then  the  derivative  of  the  functional  f in 
(4.11)  with  respect  to  t is  given  by 

f = f!  (A(-),A)  • A + 6f(A( -),A|A(X'))  . (4.17) 

/v 

• • 

Moreover,  the  functional  6f ( A( • ) , A| A(X ' ) ) , which  is  linear  in  A(X')  admits  the 
representation 

it  = f #(X,X')  • A(X')(p7p') dV' 

4(X,X ' ) • A(X ' )dv ' , (4.18) 

^ rsj  rw  m 

where  $(X,X')  is  specified  by 

*(X,X')  = * ( A(X* ) , A(X ' ) , A(X ) ,X* ,X ' ,X ) (4.19) 

/v  /v  ^ 'w  ^ rv  rv  /v 

x* 

~ 0 
A 

and  i is  a functional  of  X"  over  ft  . 

— ~ o 

Returning  to  our  main  objective  concerning  the  constitutive  equations  in 

nonlocal  elasticity,  we  recall  that  in  obtaining  the  local  field  equations  we 

have  already  tacitly  introduced  the  usual  continuity  and  differentiability 

assumptions,  which  must  be  satisfied  by  the  various  fields  and  functions. 

Moreover,  as  noted  above,  it  has  been  also  assumed  that  all  quantities  in 

(4.4)  to  (4.7)  are  continuous  functions  of  X throughout  ft  for  all  times.  We 

now  introduce  the  constitutive  assumptions  that  T),?,f ,T,p,u,X  are  continuous 

functionals  of  the  type  (4.11),  and  that  f,  in  addition,  is  a continuously 

differentiable  functional  in  the  sense  defined  above.  We  also  assume  that  p, 

• • 

is  a continuous  functional  of  A(*)  and  a continuous  function  of  A(X).  Then,  if 
t is  the  quadruplet 


19- 


$ = (m,m,n,K)  , 


/v 


we  may  write 


: = Mq  + M.  • +li.-  +M.p 
v ae  ag£  ax  * a£  ~ 


f (m0' +m  . g' +n  • x' +K  • F')dv'  . 


The  last  expression  can  be  rewritten  in  the  following  form: 

Pt  = - pTie+  l • grad  0+  (p/pQ)£  * £ 

A A 

-g8A-L-pf»v+P(i  , 


where 


and 


pT^  = p aa  + j P ,m  ,(iv ' » A = P + !*  P ,m 

ae  JR  ~ “r  ~ 


'dv' 


(p/p0}i = p If  + L p'*'dv'  » • p£=  p * + L p’"'dv’  » 


~ rt 


^ I®  • • • • 

pp,  = ! [pme'  - p 'm'0  + pm  • (grad  0) ' - p'm'  • grad  e 
R ~ ~ 

T T 

- pm  • L ' g'  + p'm'.Lg 


+ pn  • v ' - p 'n  ' • v 


+ pK  . F'  - p'K'  • F]dv ' . 

/v  m /%/ 


We  observe  that  ^ satisfies  the  condition 

A 


r A 

| p^dv  = o 


identically. 


Since  the  expression  (4.22)  for  pf  must  be  invariant  under  superposed 
rigid  body  motions,  it  follows  from  (4.22)  and  (4.24)  that  pf  has  the 


(4.20) 


(4.21) 


(4.22) 


(4.23) 


(4.24) 


(4.25) 


alternative  form 
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A |» 

pf  = (pn  - p 'n  ' )dv  ' 


(4.26) 


. A A . A A? 
p(x»f  - f ®x)  + T - T 


+ J [pn®x'  - p'n' ®x  + p KF'T  - p'k'FT-  pg ' ®m  + p'g®m']dv' 

R 

-Jj PV®i+PKF,T_-  dv'  = 0 , 


where 


A A 'p 

T - (p/p  )P  F1  - g Si 

“ ' (J  ^ /w  /v 


(4.27) 


(4.28) 


From  (4.26)  and  (4.27)  we  see  that 


r A 

I Pfdv  = 0 


°ft  ~ 


l»  A A A Am 

[p(x®f  - f Sx)  + T - T ]dv  = 0 . 

J _ ^ /W  /V 

K 


(4.29) 


(4.30) 


The  results  (4.26)  and  (4.27)  can  also  be  obtained  by  making  explicit  the  condition 
that  if  is  an  invariant  function  of  (4.4)  and  a functional  of  (4.5)  under  superposed 
rigid  body  motions,  and  by  using  the  results  (4.23)0  , , and  (4.28)  (see,  e.g.  [4]). 
For  convenience,  we  now  define  T),f,p),T  by  the  equations 


A _ a 

71  = 11  + 11  , f=f  + f , 


A - A _ 

U = n + p.  , T = T + T . 


(4.3D 


The  functions  71,f,T  satisfy  the  same  type  of  constitutive  equation  as  Tl,f,T 
and,  in  view  of  (2.16),  (2.28),  (4.29)  and  (4.30), 

f pfdv  = 0 , 


(4.32) 


f [p(x®f  - f®x)  tT-^ldv  = 0 . (4.33) 

K 

The  function  ^ has  the  same  type  of  constitutive  equation  as  p,  and  because  of 
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(2.30)  and  (4.25), 


(•  _ 

I Pudv  = 0 


Since  the  energy  equation  (3-5)  is  regarded  as  an  identity  for  all  processes 
and  since  £ and  p are  independent  of  the  rates  (4.6)  and  (4.7)  and  p,  is 
linear  of  degree  one  in  these  variables,  with  the  help  of  (4.22)  and  (4-31), 
it  follows  from  (3*5)  that 


p|9  + P • grad  9=0 


- PT]9  - l • grad  9 - pf  . v + T • grad  v + pp  = 0 


From  (4.34)  and  (4.36)  we  have 

i [ (p7]  - di v £)Q  + (pf  + div  T)  • v Jdv  - f (Tn  • v - Qi  • n)da  = 0 
"ft  ~ ~ JdR 

for  all  processes.  For  given  values  of  the  variables  (4.4)  and  (4.5)  we  may 

choose  9 and  v arbitrarily  and  independently  throughout  RUdR,  so  that 


pT)  - div  l = 0 , p.f^  + div  T = 0 in  R , 


Tn  = 0 , jt-n  = 0 on  dR 


These  are  restrictions  on  constitutive  equations  for  71, Z,T  to  be  satisfied  for 
all  processes.  Also,  from  (4.3b)  and  (4.38)  we  have 


- div(jl0-  T^v)  + pp,  = 0 


Again,  this  is  a restriction  on  the  constitutive  equations  for  Z,T,u  for  all 

^ ^ /w 


processes . 


For  the  class  of  constitutive  equations  considered  here  it  is  difficult 
to  make  explicit  the  restrictions  implied  by  (4.38)  to  (4.40).  Given  these 


restrictions  then  the  subsidiary  conditions  (4.32)  to  (4.34)  are  satisfied. 

If  the  Helmholtz  function  | does  not  depend  explicitly  on  g,g',  then 
from  (4.23)2*  (4.38)  and  (4.40)  we  see  that 

T\  = 0 , £ = 0 , pp,+  div(TTv)  = 0 , 


while  (4.28)  becomes 


A , . .A  Af 

T = (p/pn)P  F 


In  view  of  (4.38)  and  (4.39)  we  see  that  the  internal  body  force  f and  stress  T 

/v 

make  no  contribution  to  field  equations  and  yield  no  applied  stress  on  dft.  It 
may  happen  that  from  the  identities  (4.38),  (4.42)  and  (4.4l)  it  follows  that 


T = 0 , f = 0 , u = 0 


but  we  have  not  been  able  to  prove  this.  In  view  of  the  properties  associated 
with  f ,T  it  seems  natural  to  place  special  emphasis  on  the  situation  in  which 

/w  rw 

(4.43)  hold.  Given  (4.4l)  to  (4.43),  it  follows  from  (4.31)  that 

A A A A „ , .A  t /, 

TJ  = Tj  , f=f  , p = p,  , T = T = (p/p  )P  F , (i 

/% j r+s  rw  rw  U /n*  rw 


A A 


where  T],...,P  are  given  by  (4.23),  ^ 1 > (4.24)  and  (4.26).  The  field  equation 

~ ->-,3,4 

(2.14)  reduces  to 


A A 

div  T + p (b  + f ) = pv 


while,  with  the  help  of  (4.35),  equation  (3-2)  becomes 


p9T]  = pr  - div  q 


Constitutive  equations  for  T,A.  or  equivalently  T,r,  must  satisfy  the  identity 
(2.24). 


The  above  discussion  of  constitutive  equations  has  been  carried  out  without 

the  help  of  any  form  of  second  law  of  thermodynamics . When  the  free  energy 

/ * 

response  function  does  not  depend  explicitly  on  g and  g , then  *=  * and,  by 
(*+*23)2,  (1+-3l)1  and  (4.41)^,  Tl=7]  also.  Moreover,  from  an  examination  of 
(3-12),  (4.35)  and  (2.30),  it  can  be  seen  that  the  inequality  (3»ll)  is  identically 
satisfied  and  yields  no  new  information.  We  do  not  consider  here  the  case  in 
which  i|r  does  depend  on  g and  g'  and  set  aside  a discussion  of  the  conditions 
to  be  imposed  on  the  heat  flux  vector. 

The  constitutive  equation  for  i|r  and  other  results  given  in  [4]  may  be 
regarded  as  special  cases  of  those  obtained  here.  The  present  results, 
however,  differ  in  some  respects  from  those  given  in  [6,7]. 
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20.  ABSTRACT  (Continue  on  rorotoo  oldo  II  nac oooory  and  Identity  4 y 4lac4  numOot) 

^It  is  shown  in  this  pappr  that  the  conservation  laws  of  nonlocal 
continuum  mechanics  are  equivalent  to,  and  can  be  derived  from,  the 
conservation  of  energy  and  the  invariance  conditions  under  super- 
posed rigid  body  motions.  Also,  the  theory  of  nonlocal  thermo- 
elasticity is  reconsidered  in  the  light  of  recent  developments  in 
thermodynamics,  taking  invariance  conditions  fully  into  account.^ 
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